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Abstract 

Computing heavy quark- ant iquark potentials within the AdS/CFT correspondence of- 
ten leads to behaviors that differ from what one expects on general physical grounds 
and field-theory considerations. To isolate the configurations of physical interest, it is 
of utmost importance to examine the stability of the string solutions dual to the flux 
tubes between the quark and antiquark. Here, we formulate and prove several general 
statements concerning the perturbative stability of such string solutions, relevant for 
static quark-antiquark pairs in a general class of backgrounds, and we apply the results 
to M = 4 SYM at finite temperature and at generic points of the Coulomb branch. In all 
cases, the problematic regions are found to be unstable and hence physically irrelevant. 
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1 Introduction 



In the framework of the AdS/CFT correspondence [1], the computation of the static 
Wilson-loop quark- ant iquark potential in M = 4 super Yang-Mills at large 't Hooft cou- 
pling is mapped to the classical problem of minimizing the action for a string connecting 
the quark and antiquark on the boundary of AdSs and extending into the radial direc- 
tion. This approach was first employed in [2] for the conformal case and further extended 
in [3, 4] to non-conformal cases. Moreover, motivated by the recent interest in applying 
AdS/CFT ideas and techniques to calculations relevant to moving mesons in thermal 
plasmas, the finite-temperature calculations of this type have been extended in various 
ways [5]- [9]. 

Although in the conformal case these calculations reproduce the expected Coulomb be- 
havior of the potential, in the non-conformal extensions there are cases where the behavior 
of the potential is quite different from what one anticipates based on expectations from 
the field-theory side. For instance, for finite-temperature M = 4 SYM, one encounters 
multiple branches of the quark- ant iquark potential as well as a behavior resembling a 
phase transition, while for the Coulomb branch of M = 4 SYM one encounters a linear 
confining behavior at certain regions of the moduli space. To check whether these types 
of behavior are actually predicted by the gauge/gravity correspondence, one must sub- 
ject the corresponding string configurations to certain consistency checks, one of which 
is their stability under small fluctuations. 

For the conformal case, the equations of motion for the fluctuations have been first ob- 
tained in [10, 11] and result in a positive spectrum as expected. For the finite-temperature 
case, where two branches of the solution occur, the question of stability was investigated 
recently in [12] for the case of a moving quark- ant iquark pair. For this problem, the 
small-fluctuation analysis yielded a system of coupled differential equations, and numer- 
ical investigations suggested that the long string corresponding to the upper branch of 
the potential is unstable against small perturbations in the longitudinal directions. 

Here we will conduct a stability analysis, restricting ourselves to static configurations 
and pursuing an analytic treatment all the way. We will investigate perturbations in all 
coordinates including the angular ones and we will consider M = 4 SYM both at finite 
temperature and at the Coulomb branch. In this setting, the equations for the longi- 
tudinal, transverse and angular perturbations decouple and the problem is amenable to 
analytic methods. In particular, the strategy that we will employ is to write the equations 
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of motion satisfied by the fluctuations in both Sturm-Liouville and Schrodinger forms and 
to deduce the frequency spectrum of the fluctuations either from the Sturm-Liouville zero 
modes or from the Schrodinger potentials, using both exact and approximate methods. 

This article is organized as follows. In section 2, we describe the calculation of Wilson 
loops in the AdS/CFT framework via classical string configurations in D3-brane back- 
grounds and we discuss the general features of the resulting quark- ant iquark potentials. 
In section 3, we consider small fluctuations about the classical solutions, we turn the 
equations of motion of these fluctuations into one-dimensional Schrodinger problems, 
and we establish several general results for the exact and approximate calculation of the 
corresponding energy eigenvalues. In section 4, we specialize the general calculations of 
section 2 to the cases of non-extremal and multicenter D3-branes and we identify certain 
regions of the parameter space for which the potentials fail to capture the behavior ex- 
pected on physical grounds. In section 5, we apply the general stability analysis to these 
configurations and we prove that all regions where problematic behavior arises are per- 
turbatively unstable. In section 6, we summarize and conclude. Finally, in the appendix 
we discuss an interesting analogy of our problem with the standard classical-mechanical 
problem of determining the shapes of a soap film stretching between two rings and their 
stability. 

2 The classical solutions 

Our general setup refers to the AdS/CFT calculation of the static potential of a heavy 
quark- ant iquark pair according to the well-known recipe of [2] for the conformal case 
and its extensions [3, 4] beyond conformality. On the gauge-theory side, this potential 
is extracted from the expectation value of a rectangular Wilson loop with one temporal 
and one spatial side. On the gravity side, the Wilson loop expectation value is calculated 
by extremizing the Nambu-Goto action for a fundamental U-shaped string propagating 
into the dual supergravity background, whose endpoints are constrained to lie on the 
two temporal sides of the Wilson loop. Below, we give a brief review of this procedure, 
following [4], and we outline the main qualitative features of the resulting potentials. 

We consider a general metric of the form 

ds 2 = G tt dt 2 + Gy y dy 2 + G uu du 2 + G xx dx 2 + G es d6 2 + ... , (2.1) 
noting that we will use Lorentzian signature throughout the paper. Here, y denotes the 



3 



(cyclic) coordinate along which the spatial side of the Wilson loop extends, u denotes the 
radial direction playing the role of an energy scale in the dual gauge theory and extending 
from the UV at u — > oo down to the IR at some minimum value w min determined by the 
geometry, x stands for a generic cyclic coordinate, stands for a generic coordinate on 
which the metric components may depend in a particular way to be specified shortly and 
the omitted terms involve coordinates that fall into one of the two latter classes. Note 
that the metric (2.1) is diagonal and that we do not consider off-diagonal terms in the 
present paper. For the analysis that follows, it is convenient to introduce the functions 

g(u, 9) = -G tt G uu , f y (u, 9) = -G tt G yy , (2.2) 

while for the stability analysis, we also introduce 

f x (u, 9) = -G tt G xx , f e {u, 9) = -G tt G ee , h(u, 9) = G yy G uu . (2.3) 

It is useful to mention the behavior of the above functions in the conformal limit, where 
the metric reduces to that of AdS 5 x S 5 . Using the leading order expressions 

u 2 R 2 

-G tt ^ G xx ~ Gyy ~ — , G uu — — , Gee — R 2 , as u -> oo , (2.4) 

R z u z 

we see that 

g ~ h ~ 1 , f x ~ f y ~u 4 , fe^u 2 , as u -> oo . (2.5) 

In the framework of the AdS/CFT correspondence, the interaction potential energy of 
the quark- ant iquark pair is given by 

e~ iET = (W(C)} = e iS ^ , (2.6) 

where 

S[C] = J dr da ^/-det g a(3 , g a/3 = G ^d a x^dpx v , (2.7) 

is the Nambu-Goto action for a string propagating in the dual supergravity background 
whose endpoints trace the contour C. To proceed, we fix reparametrization invariance 
by choosing 

t = r , u = a , (2.8) 
we assume translational invariance along t, and we consider the embedding 

y = y(u) , x — , 9 = 9 = const. , rest = const. , (2.9) 
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supplemented by the boundary condition 

U (=4) = °° ' (2 ' 10) 

appropriate for a quark placed at y — —L/2 and an antiquark placed at y — L/2. In 
the ansatz (2.9), the constant value 6> of the non-cyclic coordinate 9 must be consistent 
with the corresponding equation of motion. As we shall see later on, this requires that 

d e g(u,e)\e=e = d e f y (u,9)\ e= 9 = , (2.11) 

which is definitely satisfied if the metric components obey identical to the above vanishing 
relations. We note that there are occasions where the latter stronger condition is not 
satisfied but (2.11) is [13]. For the ansatz given above, the Nambu-Goto action reads 

S = -£ J du^/g(u) + f y (u)y' 2 , (2.12) 

where T denotes the temporal extent of the Wilson loop, the prime denotes a derivative 
with respect to u while g(u) = g(u,9 ) and f y (u) = f y (u,9 ) are the functions in (2.2) 
evaluated at the chosen constant value 9 of 9. To seek a classical solution y c \(u), we 
proceed by reducing the problem to quadrature. Independence of the Lagrangian from 
y implies that the associated momentum is conserved, leading to the equation 

V9 + fyVc\ h 
where uq is the value of u at the turning point, f y o = f y (uo), y c \ is the classical solution 
with the two signs corresponding to the two symmetric branches around the turning 
point, and F stands for the shorthand 

F = . (2.14) 

Jy ~ fyo 

Integrating (2.13), we express the separation length as 



L = 2fJ 2 I du^- . (2.15) 

Ju Jy 

Finally, inserting the solution for y' cl into (2.12) and subtracting the divergent self-energy 
contribution of disconnected worldsheets, we write the interaction energy as 

E = - duVF-- du^fg . (2.16) 

71 Ju () 71 7u min 

Ideally, one would like to evaluate the integrals (2.15) and (2.16) exactly, solve (2.15) 
for uq and insert into (2.16) to obtain an expression for the energy E in terms of the 
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separation length L. However, in practice this cannot be done exactly, except for the 
simplest possible cases, and Eqs. (2.15) and (2.16) are to be regarded as parametric 
equations for L and E with parameter u . The resulting quark-antiquark potential E(L) 
is required to satisfy the concavity condition [14] 

§>"■ 

which is valid for any gauge theory, irrespective of its gauge group and matter content, 
and signifies that the force is always attractive and its magnitude is a never increasing 
function of the separation distance. In our case it has been shown that [4] 

dL ~2n Iy0 > U ' dl? ~ 4tt ;V 2 L'(u ) ' 1 j 

and, since in all known examples we have f Q > 0, the concavity condition restricts the 
physical range of the integration constant u to the values where the inequality 



L'(u ) < (2.19) 



is satisfied. 



In the archetypal example of [2] for J\f = 4 SYM, dual to a stack of D3-branes, Eq. 
(2.15) can be inverted for u yielding a single- valued, monotonously decreasing, function 
of L whose substitution into (2.16) gives the expected Coulomb potential, albeit with a 



coefficient which is proportional to y^JiV rather than gy M N, presumably as a result 
of strong-coupling physics. This Coulombic behavior persists in the UV limit for all 
cases involving supergravity solutions that become asymptotically AdS 5 x S 5 as u — > oo. 
However, towards the IR, the behavior of the potential is quite different and depends 
on the details of the supergravity solution. In particular, for Af = 4 SYM at finite 
temperature, dual to a stack of non-extremal D3-branes [3], as well as for the Coulomb 
branch of M = 4 SYM, dual to multicenter distributions of D3-branes [4] , one encounters 
three types of behavior, described below. 

• A solution for uq = Uq{L) exists only for L below a maximal value L c = L(uq c ) and 
is a double-valued function of L. This signifies the existence of two classical string 
configurations satisfying the same boundary conditions, called a "short" string for 
uq > u 0c and a "long" string for uq < uq c , and satisfying L'(u ) < and L'(uq) > 
respectively. Accordingly, E(L) is a double- valued function with the lower and 
upper branches corresponding to the short and long strings respectively. Although 



6 



the upper branch is energetically disfavored, the classical analysis does not clarify 
whether the long string is perturbatively unstable or just metastable and hence 
physically realizable. The latter possibility would have been physically disturbing 
since it is conflict with the concavity condition (2.17) due to the fact that U(uq) > 
for this branch. Such a behavior has been encountered for Af = 4 SYM at finite 
temperature and at certain regions of the Coulomb branch. 

• A solution for uq exists only for L below a maximal value L c and it is a single- valued 
function of L. E(L) is a single-valued function describing a screened Coulombic 
potential. However, the screening length is heavily dependent on the orientation of 
the string with respect to the brane distribution. 

• A solution for uq exists for all L and E(L) is a single- valued function of L interpo- 
lating between a Coulombic and a linear confining potential. However, a confining 
behavior is unexpected for the dual gauge theories under consideration due to the 
underlying conformal structure in combination with maximal supersymmetry. 

The last two types of behavior have been encountered for some regions of the 
Coulomb branch of M = 4 SYM. Note that in these cases it is not the concavity 
condition (2.17) that has been violated. The discrepancies arise due to the facts 
that M = 4 SYM at zero temperature is not expected to be a confining theory 
and that the screening length should not be a concept that depends heavily on the 
particular trajectory we use to compute the heavy quark interaction potential. 

In all three cases, the behavior of the interaction energy as a function of the separation is 
quite puzzling and calls for a careful interpretation. For the first case under consideration, 
it is obvious that a small-fluctuation analysis will settle the issue whether the long string 
solution is perturbatively unstable. One might then be tempted to repeat the analysis 
for the two other cases, where there is no a priori indication about the stability of the 
solutions on energetic grounds. As it will turn out, such an analysis will cast as unstable 
the parametric regions that give rise to a heavily orientation-dependent screening length 
(in the second case), and a linear behavior (in the third case). It is quite impressive, in 
our opinion, that the small-fluctuation analysis suffices to resolve the puzzles in all three 
cases. 
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3 Stability analysis 



Having described the basic features of the classical string solutions of interest, we now 
turn to a stability analysis of these configurations, with the aim of isolating the regions 
that are of physical interest and for which the gauge/gravity correspondence can be 
trusted. In this section we give a general description of the small-fluctuation analysis for 
the configurations discussed above with metric of the form (2.1) and we establish a series 
of results which will ultimately allow us to identify the stable and unstable regions by a 
combination of exact and approximate analytic methods. 



3.1 Small fluctuations 

To investigate the stability of the string configurations of interest, we consider small 
fluctuations about the classical solutions discussed above. In particular, we will be in- 
terested in three types of fluctuations, namely (i) "transverse" fluctuations, referring to 
cyclic coordinates transverse to the quark- ant iquark axis such as x, (ii) "longitudinal" 
fluctuations, referring to the cyclic coordinate y along the quark-antiquark axis, and (iii) 
"angular" fluctuations, referring to the special non-cyclic coordinate 9. To parametrize 
the fluctuations about the equilibrium configuration, we perturb the embedding accord- 
ing to 

x = 6x(t,u), y = y c i(u) + Sy(t,u) , 9 = 9 + S9(t,u) . (3.1) 

keeping the gauge choice (2.8) unperturbed by using worldsheet reparametrization in- 
variance. 1 We then calculate the Nambu-Goto action for this ansatz and we expand it 
in powers of the fluctuations. The resulting expansion is written as 

S = So + Si + S 2 + . . . , (3.2) 

where the subscripts in the various terms correspond to the respective powers of the 
fluctuations. The zeroth-order term gives just the classical action. The first-order con- 
tribution is easily seen to be 



Si = — 
2tt 



X - j dtdu 5y' + (^y^dgg + ^°^ 2 dof^j 56 



(3.3) 



1 Alternatively, we could perturb the choice of gauge by setting u = a + Su(t, a), keeping y fixed to 
t/ c i(c) as in [12]. The advantage of our choice is that the differential equation for Sy is simpler than that 
for 8u. However, as we shall see in subsection 3.2, in determining the boundary conditions the discussion 
is facilitated by using 5u instead of Sy. For details on the various possible gauges, see [15]. 
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where d$g and d$f y stand for the ^-derivatives of the full functions g(u,9) and f y (u,9), 
evaluated at 9 = 9q. The first term is a surface contribution which may cancelled out 
by adding to the classical action the boundary term which will not affect the classical 
equations of motion. In the second term, the coefficient of 59 is just the equation of 
motion of 9 and the requirement that it vanish leads indeed to the condition (2.11). Using 
this condition and performing straightforward manipulations, we write the second-order 
contribution as 



S 9 = 



2tt/ 



dtdu 



fx 



2FV2 



5x 



12 



hfxF^\„ 



9 -h. <w 2 

2F 3 / 2 y 

f 



2gf y 

h 



-5x 2 



2F 1 / 2 



Sy 2 



(3.4) 



2j pl/2 



§g/2 _ hf 9 F^ 2 5 Q 2 + 



1 



("4FV2 



q2 , fyQ Fl/2 x2 f 



69 2 



where all functions and their ^-derivatives are again evaluated at 9 = 9 . We observe that, 
by virtue of (2.11), the various types of fluctuations decouple, which greatly facilitates 
the analysis. Writing down the equations of motion for this action, using independence 
of the various functions from t, and introducing an e~ luJt time dependence by setting 



6x»(t,u) = 5x»(u)e 



-iu>t 



(3.5) 



we obtain the following linearized equations for the three types of fluctuations 



-1 


( fx d 


du 1 


yF l ' 2 du 


-1 


( 9fy d 


du 1 


\FV 2 du 




f fe d 


-1 

du ' 


V F 1 / 2 du 



gfy 



5x = , 



h 



+ u 



F 1 / 2 

2 hf s F^ 2 



5y = 0, 



(3.6) 



dig 



fyoF 1/2 2 f 



59 = . 



gf y 2FV2 - 2/2 

Therefore, the problem of determining the stability of the string configurations of interest 
has reduced to a standard eigenvalue problem for the differential operators referring to 
the three types of fluctuation. More precisely, we are dealing with differential equations 
of the general Sturm-Liouville type 

d 



d_ 

du 



p(u;u ) 



du 



- r(u;u ) > $(u) = u q(u;u )§(u) , u min < u < u < oo , (3.7) 



where the functions p(u;uo), q(u;uo) and r(u;uo) are read off from (3.6) and depend 
parametrically on uq through the function F in (2.14). Our problem then is to determine 
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the range of values of u$ for which uj 2 is negative, signifying an instability of the classical 
solution. 

Although in many cases the Sturm-Liouville description given above is sufficient for our 
purposes, in other cases it will be convenient to transform our problem into a Schrodinger 
one. To do so, we employ the changes of variables 



x 



rv/i = r 

Ju V P Ju 



dui 



h 



fy - fyO 



(3.8) 



where we note that the expression for x is valid for all three types of fluctuations. Then, 
Eq. (3.7) is transformed to a standard Schrodinger equation 

~ d 2 



dx 2 



+ V(x;u ) 



with the potential 



r p 1 / 4 d 

V = h 

q g 3 / 4 du 



1/2 



d_ 

du 



(pq) 



1/4 



^(x) = uj 2 ^(x) , 



(3.9) 



(3.10) 



expressed as a function V(u; u ) of u in the first relation and as a function V(x; u ) of x 
in the second one. The Schrodinger problem is defined in the range x G [0, x ] where x 
is given by 

/ 



x 



du. 



h 



~ fyO 



(3.11) 



and is finite, as can be verified using the asymptotic expressions (2.5) and the fact that 
p/q ~ u — u as u — > uq. The precise behavior of xo for finite values of the parameter 
u depends on the details of the underlying supergravity solution. However, using the 
asymptotic expressions (2.5) we deduce that 

r(i/4) 2 



X(U ) 



as -Uq — > oo 



(3.12) 



4V27T Uq 

The fact that the size of the interval is finite implies that the fluctuation spectrum is 
quantized and the fact that it becomes narrower in the UV implies that there cannot be 
any instability in the conformal limit of the theory. With the transformation (3.8) the 
UV and IR regions are mapped to the regions near x = and near x = xo respectively. 
In the Schrodinger description, our problem is to determine the range of uq for which the 
ground state of the potential in (3.10) has negative energy. Note that, in general, the 
first of (3.8) does not lead to a closed expression for u in terms of x and therefore it is not 
always possible to write down the potential as an explicit function of x. Nevertheless, in 
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many cases the expression for the potential as a function of u provides useful information 
by itself (for example, if V(u;uq) is manifestly positive then the energy eigenvalues uo 2 
are always positive and no instability occurs), while its transcription in terms of x may 
be easily accomplished when performing perturbation-theory calculations. 



3.2 Boundary conditions 

To fully specify our eigenvalue problem, we must impose appropriate boundary conditions 
on the fluctuations at the UV limit u — > oo (x = 0) and at the IR limit u — u (x — x ). 
Starting from the Sturm-Liouville description, the boundary condition at the UV is 
particularly easy to determine by looking at the singularity structure of the differential 
equations in (3.6) as u — > oo. Indeed, at this limit we find 

$" + 1$' = 0(u~ 4 ) , as u — > oo , (3.13) 

u 

which implies that u = oo is a regular singular point and that the two independent 
solutions of any equation in (3.6) have the form 



?1=0 



1 oo 

d\nu + — V d n u~ n . (3.14) 



u 

n=0 



For the first solution, which is regular at infinity, we use the freedom to assign to the 
constant Co any value we want to set it to zero. For the second solution, which blows up 
at infinity, we choose d — 0. Hence we impose the boundary condition 

$(«) = , as u — > oo . (3.15) 

The nature of the boundary condition at the IR needs special attention due to the fact 
that we have to glue the fluctuations around the upper and lower branches of the classical 
solution y c i, corresponding to the two different signs in (2.13). The point u = u is again 
a regular singular point and the two independent solutions admit expansions of the form 2 



$! = (u-u y 



CO 



n=l 



$2 = d + J2dn(u-u ) n , (3.16) 



2 We assume that u > u m i n strictly. In case of equality the singularity structure is different and will 
be examined in general and in the examples of particular interest below. 
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where p = 1/2 for the fluctuations along x and 9 and p = —1/2 for those along y. We 
demand that, on the two sides of the classical solution, the fluctuations and their first 
derivatives with respect to the classical solution be equal. However, this should be done 
not using the classical solution y ch but that combined with its perturbation as in (3.1). 
Equivalently, defining a new coordinate u as 

u — u + 5u(t, u) , 5u(t,u) — — ^yjt ,v ) ; (3.17) 

y c i( u ) 

the classical solution is not perturbed at all. This redefinition does not affect the x- 
and ^-fluctuations since they have trivial classical support and we keep only linear in 
fluctuations terms. However, since y' cl ~ (u — -u ) _1 ^ 2 near u — M o? the fluctuation 5u 
has an expansion as in (3.16), but with p = 1/2. Since, around the matching point, the 
two branches of the classical solution differ by a sign while the fluctuations are equal, we 
have 

au U=UQ 

where <3> refers to the Sx, 89 and 8u fluctuations. Recasting everything in terms of the 
original fluctuations we see that, in the expansions (3.16), the coefficient c = for the 
x and 9 fluctuations and do = for the y ones. Equivalently, 



8x,80\ (u-u ) 1/2 & = , as u 



u 



i 



5y: $ + 2{u - uq)& = , as u -> m+ . (3.19) 

These boundary conditions for the Sturm-Liouville function $ can be transcribed to 
boundary conditions for the Schrodinger wavefunction \l> by using the asymptotic rela- 
tions u ~ 1/x for u — > oo and u — u ~ (x — x) 2 for -u — > u and x — > Xo, which 
follow from (3.8) and the asymptotic behavior of the p- and g-functions. After some 
manipulations, we find that for all types of fluctuations we must have 

(0) = , *'(xo) = , (3.20) 

which correspond to a Dirichlet and a Neumann boundary condition, in the UV and IR 
respectively. 



3.3 Zero modes 

The method described above would in principle allow us to determine the regions of stabil- 
ity of the solutions, provided that the relevant Sturm-Liouville or Schrodinger problems 
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could be solved exactly. However, in many cases, these problems are quite complicated 
and the spectrum is impossible to determine. On the other hand, it turns out that we can 
obtain useful information by studying a simpler problem, namely the zero-mode problem 
of the associated differential operators. In what follows, we prove that (i) transverse zero 
modes do not exist, (ii) longitudinal zero modes are in one-to-one correspondence with 
the critical points of the function L(uq), and (iii) the angular zero-mode spectrum can 
be obtained to good accuracy by approximating the corresponding Schrodinger potential 
by an infinite square well. These results are crucial for the stability analysis of section 5. 

3.3.1 Transverse zero modes 

We consider first the case of the transverse fluctuations. The zero mode solution obeying 
(3.15) is up to a multiplicative constant 

du / gf y 

fxyfy — fyO 

= ^v^-sf^VR^M^ (3 - 21) 

Jxjy Ju \ Jxjy I 

Ju V Jxjy J JxO y JyO 

where in the second step we performed a partial integration. The zero mode solution 
exists provided that we can satisfy the boundary condition (3.19), i.e. make the coefficient 
of (u — Wo) 1 / 2 vanish for some value of the parameter u . It is easy to see that this is not 
possible and so transverse zero modes do not exist. Therefore, if the lowest eigenvalue 
of the Schrodinger operator corresponding to the transverse fluctuations is positive for 
some value of u , it will stay positive throughout. Since, as it will turn out, this is indeed 
the case for all situations under consideration, the classical solutions are stable under 
transverse perturbations. 

3.3.2 Longitudinal zero modes 

We next turn to the longitudinal fluctuations, for which we will prove the powerful result 
that normalizable zero modes exist only at the values of uq where the length function 
L(uq) in (2.15) has a critical point, i.e. L'(uq) = 0. The importance of this result lies in 
the fact that it implies that a point where the lowest eigenvalue u 2 changes sign must 
necessarily be a critical point of the length function. 
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To establish the connection between longitudinal zero modes and critical points of the 
length function, we begin by writing the longitudinal zero-mode solution obeying (3.13). 
Up to a multiplicative constant, we have 

\/!Pfy 



Sy 



J 

J u 



Jy V •'.</ JyO 



du 



a/ fy ~ fy 

9ofy ° ( U -u )- i / 2 +2 r 

J Ui 



-.d„ 



du 



f 
J y 



d u 



(3.22) 



fyO Ju y/fy TyO \ fy 

Then, the boundary condition (3.19) implies that this mode exists only if 

du „ / \fgfy 



Juo 







(3.23) 



\f fy - fyO V f'v 

We next differentiate the length function (2.15) with respect to u . The result is easily 
seen to be 



L'(u ) 



;y0 



Jun 



fyO Juo 
fyO 



du 



Wy 



f 



(fy ~ fyo) 3/2 
du 



_ 9f l/2 Fl/2 
Z JyO f 

Jy 



U=UQ 



f, 



fyO Ju 

Tn lim 

u—>uo 

yo 



a/ fy - fyO 



d u 



f 
J y 



(3.24) 



jpl/2 ( dug fyO _ fyO 
dufy fy 



where as before we have performed a partial integration. The last line is zero, so that 



L'{uq) 



fy 



yO 



du 



d u 



^y 



(3.25) 



^ fyO Ju \f fy ~ fyO \ fy 
which is manifestly finite, in contrast with the first line in (3.24) where both terms 
diverge. Comparing (3.23) and (3.25) we see that at the extrema of L(u ) there is a zero 
mode, as advertised. As a result, given the critical points of the length function, we can 
determine the sign of the lowest eigenvalue u 2 in all regions by expanding the Schrodinger 
potential in uq about each critical point uq c and determining whether uj 2 changes sign 
there. 

There is an alternative way to understand the occurrence of the critical value u 0c in the 
above analysis. It is not difficult to prove that the Schrodinger potential has the following 
behavior at the extreme values of u, namely that 



V(u; u ) = 2u 2 , as u — > oo 



(3.26) 
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and 

%W= ^) = & + ^-|^g + g. (3,7, 

Since the potential rises from a minimum value to infinity in the finite interval x G [0, x ], 
where x is given by (3.11), the spectrum of fluctuations is discrete. In addition, if its 
minimum value V is negative enough the potential can very well support bound states 
with negative energy and this happens whenever (3.23) has a solution for some values of 
the parameter uq. The largest of these values is woe in which the lowest energy eigenvalue 
becomes zero. 



3.3.3 Angular zero modes 

We finally consider the angular zero modes which, unlike the two other types of zero 
modes, cannot be explicitly written in an integral form due to the presence of the "restor- 
ing force" term in the corresponding Sturm-Liouville equation at the third line of (3.6). 
On the other hand, in the Schrodinger description, we can obtain the full fluctuation 
spectrum to quite high accuracy by approximating the corresponding potential by an 
infinite well. The procedure is described below. 

We first note that the behavior of the angular Schrodinger potentials at the limits u — > oo 
and u = uo in all our examples is as follows 

1 



and 



Vooiuo) = V(oo; u ) = - lim u 2 d 2 e g , (3.28) 

V (u ) = V(u ; u ) = l 9 -^ d U0 (!*&) + l-fir%f» ■ (3-29) 

8 n oJeo \9oJyo/ ^'T'Ojeo 

where the limit in the first equation is finite. Since the potential rises from a minimum 
value to infinity in the finite interval x G [0, x ], where x is given by (3.11), we may 
approximate it by an infinite well given by 

V; pP ,«4^ + y (3.30) 
I oo , otherwise J 

With the boundary conditions (3.20) the energy spectrum is given by 



vr 2 (2n + l) 2 1 



u n ( u o) = 4 ^ o ( Mo )]2 + 2 [^(^p) + V o°( u o)\ , n = 0, 1,... . (3.31) 

Note that we have taken the average value of the values of the potential at the two 
extreme values of u. In order for a solution to exist, the above average must be negative 
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at least in a finite range of values for uq. The critical value uq c is then obtained by 
numerically solving the equation Uq(uq) = 0. We note that different choices for V^ pprox 
are possible (e.g. the average value of the potential itself) but, conceptually, they do not 
offer something new. 

The advantage of the method just described is that it gives approximately the entire 
spectrum of fluctuations and its dependence on the parameter u without the need for 
more sophisticated numerics. In particular, we may plot the lowest eigenvalue oJq as 
a function of uq and easily check whether it is monotonously increasing function, as it 
is in all of our examples. In addition we will see that, in one case, the infinite- well 
approximation is in fact exact. Finally, we will show that there is a hierarchy of different 
critical values Woc,n, n = 0, 1, ... , with u 0c = -u c,o an d w 0Ci „ > u oc,n+i, at which the 
corresponding eigenvalues become zero and below which they become negative. 

At this point one might wonder whether this infinite-well approximation could be applied 
to the longitudinal fluctuations, using, for example, (3.27) for the value of the potential at 
the bottom of the well, to yield in addition the full approximate spectrum of fluctuations. 
For our specific examples, it can be shown that, although such a choice gives a reasonable 
approximation to the result for u 0c following from (3.23), it also results in a tower of 
negative eigenvalues as we lower u , which is in conflict with our general conclusions 
following from (3.27). The discrepancy is traced to the fact that the true potential 
varies significantly during the interval [0, x ] and thus supports fewer negative-energy 
states than the infinite- well potential used to approximate it. One may of course devise 
phenomenological potentials that mimic the expected behavior, but there is no need to 
do so in the present context. 

3.4 Perturbation theory 

As stated earlier on, the stability of the solutions against longitudinal perturbations will 
be determined by identifying the critical points of L(uq) and solving the Schrodinger 
equation for small deviations of Uo from each critical point u 0c by means of perturbation 
theory. In doing so, we must note that the parameter uq enters into the problem not 
only by appearing explicitly in the potential V(x;uo), but also by controlling the size 
of the interval Xo in which the problem is defined. In what follows, we will present the 
perturbation-theory formulas that are appropriate for such a case. 

To this end, we consider a more general Schrodinger problem defined on the interval 
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[0, xq] with a potential that depends on a parameter uq, for which we know that for some 
uq = uq c and xq = xq c the Schrodinger equation (3.9) admits a single solution with a 
given eigenvalue u 2 (equal to zero in our cases of interest). We would like to determine 
the correction to the energy eigenvalue when uq and Xq deviate from m 0c and x 0c . An 
easy computation shows that the potential can be written as 



V(x;u ) 
5V(x) 



V(x; u 0c ) + 5V(x) , 
5x 



%0c 



2V(x; u 0c ) + xd x V(x; u 0c ) + Su d UOc V(x; u 0c ) + 



(3.32) 



where x G [0, xo c ] while Suo = uq — uq c and 5xq = xq — xq c . Then, a careful computation, 
keeping track of boundary terms, gives the energy shift 



Su 2 



/ dx \^(x)\ 2 5V(x) 
Jo 



Sxo 

x 0c 



2cu 2 + 



+ w*') - x |^'| 2 + (c 2 _ v)\V 







^Oc ^ 


to 










o J 



rxoc 

+ 5u / dx \^(x)\ 2 d UOc V(x; u 0c ) , 
Jo 



(3.33) 



where the term (uj 2 — V)\^\ 2 can be alternatively written as — \ (ty'ty* + using the 

Schrodinger equation. In deriving this, we have to use the boundary conditions so that 
the unperturbed Hamiltonian is Hermitian, which is always the case when the boundary 
conditions are Dirichlet, Neumann or a linear combination thereof. 

To apply this result to our case, we set u 2 = and we note that the parameters Xo and 
uq are related by (3.11). Then, the energy shift is found to be 



5uj 2 = Suo 



x' (u 0c )V(x 0c ;u c)\^(xoc)\ 2 + / dx \^(x)\ 2 d UOc V (x; u 0c ) 



■r<)c 



(3.34) 



with xoc = xo(uoc) and where we have used the boundary conditions (3.20). Regarding 
the first term, an explicit computation, using (3.11), gives 



i>oo 

x' (u ) = fy / du 

J uo 



y/fy- fyO 



(3.35) 



Turning to the second term we must note that the Wo-derivative acts on the potential 
V(x; Uq) whose explicit form is not known. To transcribe this expression into one involv- 
ing the potential V(u; uo), for which explicit expressions are available, we must also take 
into account that when varying u while keeping x constant, u varies as well. That is to 
say, d U0 V(x; u ) is actually the "connective" derivative 



du 

d U0 V{x;u ) = d U0 V(u;u ) + ^— d u V(u;u ) , 

OUq 



(3.36) 
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where the second term can be evaluated with the aid of the expression 



Note that, in the case where (3.34) gives 5uj 2 = 0, we have to go beyond first-order 
perturbation theory to determine whether the energy eigenvalue changes sign. Luckily, 
no such behavior occurs in our examples. 



4 D3-brane backgrounds: The classical solutions 

In this section, we review the behavior of the quark-antiquark potentials emerging in 
Wilson-loop calculations for non-extremal and multicenter D3-brane backgrounds. The 
purpose of this review is to identify the three types of problematic behavior referred to at 
the end of section 2, which motivated our stability analysis. Since we will work with the 
Nambu-Goto action, we need mention in the expressions below only the metric and not 
the self-dual five-form which is the only other non-trivial field present in our backgrounds. 



4.1 Non-extremal D3-branes 

We start by considering a background describing a stack of N non-extremal D3-branes. 
The field-theory limit of the metric reads 



a 2 u2 



-[\-^)de + dxi 



' • ( dir • ,/(>r • . ( J.1) 



where the horizon is located at u — fj, and the Hawking temperature is T = -J^. This 
metric is just the direct product of AdS 5 -Schwarzschild with S 5 and it is dual to Af = 4 
SYM at finite temperature. For the calculations that follow, it is convenient to switch to 
dimensionless variables by rescaling all quantities using the parameter /i. Setting u — > fiu 
and Uq — > fiuo and introducing dimensionless length and energy parameters by 

L^—L, E^^E , (4.2) 

we see that all dependence on \x and R drops out so that we may set /i — > 1 and R — > 1 
in what follows. The functions in (2.2) and (2.3) depend only on u (reflecting the fact 
that all values of 9 are equivalent) and are given by 

g(u) = 1 , f y (u) =u 4 -l, (4.3) 
18 



(a) (b) 

Figure 1: Plots of L(u ) and E(L) for the non-extremal D3-brane. In this and subsequent 
plots, the various types of lines correspond to stable (solid dark), metastable (solid gray) 
and unstable (dashed gray) configurations, with the stability determined by the analysis 
of section 5. 

and 

/x(«) = « 4 -l, f e (u) = ^±, h (u) = -£-, (4.4) 



u 



respectively. From this it also follows that the equations of motion (2.11) for 9 are 
identically satisfied for all values of 9 . 

Let us now evaluate the quark-antiquark potential according to the guidelines of section 
2. The integrals for the length and energy are given by [3, 4] 

du _ 2 v / 2vr 3 / 2 y/u% - 1 ^ fl 3 5 1 

, uo v /( M 4_ 1)(M 4_^y- T(l/4f 

and 



! v u »- 1 / 7PrwTi = rw— ,F, i?rf<] ' (4 ' 5) 



e -i *• i - o - (--d = -fS5-^ (4 ~i' ? 4> +l ■ (4 - 6) 



where 2-^1(0, b, c; x) is the hypergeometric function and u ^ 1. For u ^> 1 the behavior 
is Coulombic while at the opposite limit, tt — * 1, we have the asymptotics 

-^(m^-f), .^(m-^-x-l). (4,) 

The function L(uq) has a single global maximum, which together with the corresponding 
maximal length and energy are given by [4] 

m c ~ 1.177 , L c ~ 0.869 , E c ~ 0.093 . (4.8) 

For L > L c , only the disconnected solution exists. For L < L c , Eq. (4.5) has two 
solutions for uq, corresponding to a short and a long string respectively and, accordingly, 
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E is a double-valued function of L. Moreover, there exists another value of the length, 
given in our case by L c ~ 0.754, above which the disconnected configuration becomes 
energetically favored and the short string becomes metastable, which can be used as an 
alternative definition of the screening length. The behavior described above is shown in 
the plots of Fig. 1 and, provided that the upper branch of E(L) is physically irrelevant, 
corresponds to a screened Coulomb potential. 



4.2 Multicenter D3-branes 

We now proceed to the case of multicenter D3-brane distributions. These were first 
constructed as the extremal limits of rotating D3-brane solutions [16, 17] in [18, 19] 
and belong to the rich class of continuous distributions of M- and string theory branes 
on higher dimensional ellipsoids [20]. These distributions have been used in several 
investigations within the AdS/CFT correspondence, starting with the works of [21, 4]. 
Here, we will concentrate on the particularly interesting cases of uniform distributions of 
D3-branes on a disc and on a three-sphere. 



4.2.1 The disc 

The field-theory limit of the metric for N D3-branes uniformly distributed over a disc of 
radius r reads 

,u 2 + rl cos 2 9 



ds * = ff-l/2(_ rft 2 + d£ 3 )+g l/2 " T'O^ du 2 

v 4 ' u + rl 

+ H 1 ' 2 [(u 2 + r 2 cos 2 9)d9 2 + r 2 cos 2 9dQ 2 3 + (u 2 + r 2 ) sin 2 9d(P 2 } , (4.9) 



where 



H u 2 (u 2 + r 2 cos 2 9) ' (410) 
while dQl is the S 3 metric. Since the only scale parameter entering into the supergravity 
solution is r it is convenient to measure lengths and energies using this as a reference 
scale. Setting u — > r u and u — > r u and introducing the dimensionless length and 
energy parameters by 

L^—L, E^^E. (4.11) 

r 7T 

all dependence on r and R drops out so that we may set r — > 1 and R — > 1 in what 
follows. The functions in (2.2) and (2.3) now depend on 9 and read 

g(u, 9) = u2 + ™* 9 , f y (u, 9) = u 2 (u 2 + cos 2 9) , (4.12) 
U z + L 
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(a) (b) 
Figure 2: Plots of L(u ) and E(L) for the disc at 9 = 0. 



and 



7/2 _j_ pQg2 Q 

f x (u, 9) = u 2 (u 2 + cos 2 9) , f e (u, 9) = u 2 + cos 2 9 , h(u, 9) = T —— , (4.13) 

respectively. The conditions (2.11) are satisfied only for 9 = and 9 = vt/2, which 
correspond to trajectories orthogonal to the disc and lying on the plane of the disc, re- 
spectively. To evaluate the quark- ant iquark potential, we examine these two trajectories 
in turn. 

• 9 = 0. For this case, the integrals for the dimensionless length and energy read [4] 

du 



L = 2u ^u 2 + l j 

Ju 



«V / (« 2 + 1 )(« 2 -«o)(" 2 + "o+ 1 ) 



*^[n(V,*)-K(*)] 





(4.14) 



and 



E 



du 



HQ 



U 2 + 1 



U ( 



- 1 



- u 



{u 2 -u 2 )(u 2 + u 2 + l) 

= \j2ul + 1 [k' 2 K(k) — E(fc)] , (4.15) 

where K(fc), E(fc) and II(ai, fc) denote the complete elliptic integrals of the first, second 
and third kind respectively and 

M 



k 



k' = , 



(4.16) 



V 2u l + 1 ' 

are the modulus and the complementary modulus. For uq >> 1, the behavior is Coulombic 
as before and at the opposite limit, uq — > 0, we have the asymptotics [4] 



L ~ 7r(l — w ) , 
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4 -o> 



(4.17) 




u 



(a) 



(b) 



Figure 3: Plots of L(u ) and E(L) for the disc at 9 = ir/2. Note that the screening 
length is less than the the one for 9 = by a factor of 2. 



Now, L{uq) is a monotonously decreasing function and hence its global maximum is at 
Moc = given by L c = ir. For L > L c , only the disconnected solution exists, while for 
L < L C) Eq. (4.14) has a single solution for u and £ is a single- valued function of L. 
This behavior is shown in the plots of Fig. 2, and it corresponds to a screened Coulomb 
potential. 

• 9 — 7r/2. Now, the integrals for the dimensionless length and energy read 

du 



2ui 



UQ 



2u\ 



and 



E 



roc 
J UQ 



duu 



u 



n 



ut 



Uq + u^ 



,k)-K(k) 



^/{u 2 + \){u A -ul) VvtTT 

: K(k) - sju~l + u\ E(fc) + 1 , 



f-UQ 
JO 



duu 



where now 



k 



u\ — u\ 
u\ + u\ 



k! = Vl^W 



(4.18) 



(4.19) 



(4.20) 



and u > (w<) denotes the larger (smaller) between u and 1. At the limit u — > 0, we have 
the asymptotics [4] 



t 71 
L ~ — 

2 



l-«n 



hw - 1 



hw - - 



(4.21) 
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The behavior is qualitatively similar to the previous case with the maximal length being 

L c = \ ■ (4.22) 

Comparing the 9 = and 9 = n/2 cases we note that, although the qualitative behavior 
of the potential is the same, the expressions for the screening length differ by a factor of 2. 
This factor is quite large as one expects that the orientation of the string configuration will 
have a mild effect on physical observables of the gauge theory. The apparent discrepancy 
will be resolved by our stability analysis. 

4.2.2 The sphere 

The field-theory limit of the metric for N D3-branes uniformly distributed over a 3-sphere 
of radius r reads 

rfs 2 = H -l/2 { _ df 2 + df 2) + jjl/2 ^ - r 2 COS 2 9 ^ 

U l — 7q 

+ H 1/2 [{u 2 - r 2 cos 2 9)d9 2 + u 2 cos 2 9dVi\ + (u 2 - r 2 ) sin 2 9d(j) 2 ] , (4.23) 

where 



u 2 (u 2 — r 2 cos 2 

Note that this solution is obtained from the disc by taking r 2 — > — Tq. Employing the 
same rescalings as before, we write the functions in (2.2) and (2.3) as 

it 2 — rns 2 9 

9^ e ) = U 2 -i > />> °) = «V - co§2 6) ( 4 - 25 ) 

and 

f x (u,9) = u 2 (u 2 -cos 2 9) , f e (u, 9) = u 2 - cos 2 9 , h(u,9) = ^3^, (4.26) 

XL X 

respectively, and the conditions (2.11) are again satisfied only for 9 = and 9 = n/2. 
We examine these two cases in turn. 

• #o = 0. For this case, the integrals for the dimensionless length and energy read [4] 

du 



2uq^uI - 1 / 



u Uy/ {u 2 - l){u 2 - ul)(u 2 +Uq-1) 

2u k' 



Un 



[U(k' 2 , k) - K(k)} (4.27) 
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(a) 



U -1.5 




(b) 



Figure 4: Plots of L(u ) and E(L) for the sphere at 9 = 0. Note the appearance of two 
branches in the plot of E(L). 



and 



E = 



du 



UQ 



u 2 - 1 



U \ 



U 2 — ul)(u 2 + Uq — 1) 



- 1 



K - 1) 



^2u 2 - 1 [k' 2 K(k) - E(A;)] + 1 , 



where 



k 



Uq 



k! = VT^¥ 



(4.28) 



(4.29) 



For uq ^> 1, the behavior is Coulombic, whereas in the opposite limit, u — > 1, we have 
the asymptotics [4] 

m - 1 



L ~ v/2( Mo - 1) 



In 



M - 1 



- 2 



In 



Uq - 1 



-3 



(4.30) 



The function L(w ) has a single global maximum. Its location, its value and the corre- 
sponding value of the energy are [4] 



u 0c ~ 1.125 , L c ~ 1.002 , £ c ~ 0.063 . 



(4.31) 



For L > L c , only the disconnected solution exists. For L < L c , Eq. (4.27) has two 
solutions for u and E is a double-valued function of L. This behavior is shown in the 
plots of Fig. 4 and, discarding the upper branch of E(L), it corresponds to a screened 
Coulomb potential. 

• 6 = 7r/2. Now, the integrals for the dimensionless length and energy read 



L = 2u\ 



du 



... Ma/(m 2 - 1)(« 4 - Mq) 

24 



V2 

Uq 



U[-,k) -K(fc) 



(4.32) 
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(a) 



(b) 



Figure 5: Plots of L(u ) and E(L) for the sphere at 9 = ir/2. Note the appearance of a 
confining potential for large values of L. 



and 



E = 



duu 



u 



'u Vu 1 

where now 



u' 1 



k 



"o 



duu 



u 



Vu^ 



-± [K (k) - 2E(fc)] . (4.33) 



2ul 



k' = VT^ 2 



For u ^> 1, the behavior is Coulombic, whereas in the opposite limit, u 
the asymptotics [4] 



(4.34) 
1, we have 



L ~ 



V2 



In 



16 



u 



- 2 v / 21n(l + V2) 



E 



2V2 



In 



16 



u 



-4 



(4.35) 



Now, L(uo) is a monotonously decreasing function which approaches infinity as u — > 1 
and zero as Mo 00 and hence no maximal length exists, Eq. (4.32) has a single solution 
for u , and E is a single-valued function of L. Therefore, it seems that no screening 
occurs but, instead, from Eq. (4.35) we have a confining potential 



£ ~ — , 

2 ' 



for L > 1 , 



(4.36) 



whose appearance is quite puzzling, given that the underlying theory is Af = 4 SYM. 

Comparing the 9 = and # = tt/2 cases, we note that the qualitative behavior of the 
potential appears to be completely different: in the first case the potential is a double- 
valued function of L which exists up to a maximal length, whereas in the second case the 
potential exists for all values of the length and interpolates between a Coulombic and 
a confining potential. One does not expect such major qualitative differences to occur 
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by a simple change of the orientation of the string and, moreover, the appearance of a 
confining potential is unexpected. Both of these issues will be resolved by our stability 
analysis. 



5 D3-brane backgrounds: Stability analysis 

In this section we apply the stability analysis developed in section 3 to the string solutions 
reviewed in section 4. We find that for the non-extremal D3-brane and the sphere with 
#o = we have a longitudinal instability corresponding to the upper branch of the energy 
curve, while for the disc with 6*0 = and for the sphere with 9 = ir/2 we have angular 
instabilities towards the IR, even though the potential has a single branch. 



5.1 The conformal case 

As a simple example, and as a consistency check, let us first consider the conformal case, 
corresponding to the \i — > or r — > limit of any of the above solutions. In this case, 
the Schrodinger potentials read 



V x (u;u ) = 2u 2 , 

V e (u;u ) = 0, (5.1) 



„4 _ „,4 

V y (u;u ) = 2 



U — Uq 



u z 

where the subscript indicates the corresponding fluctuating variable. The change of 
variables (3.8) explicitly gives 

-H^S' (6 - 2) 

and, hence, x G [0, Xq] with 

T(l/4) 2 1 , x 

x = 1 — . 5.3 
4^ u 



Since all Schrodinger equations are defined in a finite interval with positive-definite po- 
tentials, the corresponding energy eigenvalues are positive and so there is no instability. 
The equations corresponding to (3.7) for the transverse (Sx) and longitudinal (Sy) fluctu- 
ations are given by Eq. (7) of [10] and Eq. (16) of [11] respectively. The comparison for 
5x is immediate if we change variables as u — 1/z and also rename z m = l/u , where z 
and z m are the variable and parameter used in [10]. For Sy, the comparison involves the 
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same change of variables as before and 5y = u^u 2 (w 4 — -u 4 ,)" 1 / 2 5y , where Sy is the variable 
called u in [11]. In addition, since the independent variable in Eq. (16) of [11] is y (x 
in their notation) we should also use the differential relation dy = Wo/w^w 4 — -u 4 ) _1//2 rfM, 
resulting from the classical equation of motion (2.13). 



5.2 Non-extremal D3-branes 

We next proceed to the case of non-extremal D3-branes, where we recall that the potential 
energy is a double-valued function of the separation length. The Schrodinger potentials 
for the three types of fluctuations are given by 



V x (u;u ) = 2 



u 8 -u 4 



u 



6 



V (u;uo) = 0, (5.4) 

V v(^) = 2 -^-^ . 

The value of the endpoint xq is found using (3.11) and reads 

x - r(V4)2 1 F ( l 1 3 - M (55) 

x °-772T^ 2Fl U'4'i'^J ' (5 - 5) 

with the behaviors (3.12) and 

x (u ) ~ ~ ln( Mo - 1) + 7r+ g" 64 + C(m - 1) . (5.6) 

Since V x is positive for all values of uq, the solution is stable under transverse perturba- 
tions, in accordance with our general result in section 3.3. Also, Vq is identically zero, 
which means that, given that Xq is finite, the spectrum is positive-definite and the solu- 
tion is stable against angular perturbations as well. On the other hand, V y starts from a 
negative value at u = uq given by 

V y0 = a > -8^V y0 <0 (5.7) 

and behaves as in (3.26) at u — > oo. To examine the occurrence of instabilities under 
longitudinal perturbations, we apply the results of section 3.3 and the perturbation- 
theory formulas of section 3.4 to our problem. We find that there is a longitudinal zero 
mode at uq = uo c where uo c is found by numerically solving the equation (3.23) which, 
for our case, reads 

i? 2Fi G' \ i ; 4) = 2Fi G' * ^ 4) ■ (5,8) 
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This has precisely one solution given by the first of (4.8), i.e. uq c — 1.177. To find the 
change in the energy eigenvalue u 2 as we move away from uq c , we first use (3.22) to find 
the explicit expression for the longitudinal zero mode 

1.463 „ /3 3 1 7 ut 1 



Oc 



where Fi(a ; bi, b 2 , c ;x,y) is the Appell hypergeometric function and the normalization 
has been fixed from the corresponding Schrodinger wavefunction which reads 





~ U 0c 




- 1 



Then, using Eq. (3.34), we find that the change in energy eigenvalue is given by 

5cu 2 ~ (51.9 + 29.6 - 20.1) 5u = 61.4 5u , (5.11) 

where each term inside the parentheses represents the contribution of the corresponding 
term in (3.34). Hence, as we move to the right (left) of uq c , the eigenvalue u 2 becomes 
positive (negative) and therefore the long string is unstable under longitudinal perturba- 
tions. This is what is expected from the energetics of these configurations, shown in Fig. 
1(b). 

5.3 Multicenter D3-branes 

We finally turn to the case of multicenter D3-branes, where all types of problematic 
behavior discussed in section 2 appear. In what follows, we present the results of the 
stability analysis for all cases discussed in section 4.2. 

5.3.1 The disc 

For the disc distribution, we recall that both allowed orientations of the string lead to a 
screened Coulomb potential, but the screening lengths differ by a factor of 2. The results 
of the stability analysis for the two orientations are as follows. 

• 0q — 0. In this case, we have the Schrodinger potentials 

8m 8 + 18m 6 + 11m 4 - [2u 2 (u 2 + 1) - l]u 2 + u 2 (u 2 + 1) 



V x (u;u ) = 



4u 2 (u 2 + r 



2 



2u e + + [Qu 2 (u 2 + 1) - l]u 2 + 3u 2 (u 2 + 1) 
V *(u,u ) = 4u 2 (u 2 + l) 2 ' (5 ' 12) 

T . , . 8u 8 + 18m 6 - [8u 2 (u 2 + 1) - 11K - [6u 2 (u 2 + 1) - l]u 2 - 3u 2 (u 2 + 1) 
Vy[U]Uo) ~ 4 U 2( U 2 + 1)2 > 
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and the value of the endpoint x reads 

Xo = ^2k' 2 - 1 K(fc') , (5.13) 

where k' is the complementary modulus defined in (4.16). We have the general behavior 
(3.12) and 

x (uo) ~ -\n^ + O(u 2 \nu ) . (5.14) 

Since V x and V y are positive for all values of u , the solution is stable against transverse 
and longitudinal perturbations, in accordance with the general results of section 3.3 in 
the absence of a maximal length. On the other hand, Vg is negative throughout the whole 
range of uo, with its values at u — uq and u — > oo given by 

Vfco = -\ ■ (5.15) 

To examine the occurrence of instabilities, we use the infinite-well approximation of (3.30) 
and we examine the behavior of the lowest eigenvalue u 2 by plotting it as a function of 
uo- We find that it is an increasing function starting at negative values for uq = and 
changing sign at a critical value u 0c . This value, and the corresponding maximal length 
and energy, are given by 

u 0c ~ 0.48 , L c ~ 0.59vr , E c ~ -0.15 . (5.16) 

Therefore, when the separation distance of the quark-antiquark pair becomes larger than 
the value L c given above, small fluctuations in 9 destabilize the corresponding classical 
solution and the resulting potential should not be trusted. The true screening length is 
thus given by the second of (5.16) and turns out to be comparable to that for 9 = ir/2. 
Moreover, within the infinite-well approximation, we may show that for u < u 0c more 
states become negative. Indeed, setting w 2 (u ) = and using the limiting behavior 
(5.14), we find that the value -u c,n in which the n-th energy eigenvalue becomes zero is 
given by the formula 

u 0c , n ~ Ae~< 2n+1 ^ 2 , (5.17) 

valid in practice for all n ^ 1. 
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• 6*o = 7r/2. Now, the Schrodinger potentials are given by 



V x (u;u ) 



2u 6 + u 4 + u 4 
m* 



V e (u;u ) 



1 , 

2m 6 + u 4 - 2m 4 m 2 - 3u 4 
u 4 



(5.18) 



Vy(u;u ) 



and the value of the endpoint xq reads 




1 



(5.19) 



X = 



where k is the modulus defined in (4.20). Its behavior is given by (3.12) and 




(5.20) 



Since V x and V$ are manifestly positive, the solution is stable under transverse and angular 
perturbations. Also, although V y has a negative part below some value for u, the fact 
that no critical points of the length exist in this case implies that the solution is stable 
under longitudinal perturbations as well. 

The upshot of this analysis is that the screening lengths for the two different extreme 
angles for which the heavy quark potential can be computed become practically the same 
which is a requirement for the notion of a screening length to make physical sense. It 
is natural to expect that, if the system starts with O = and L c < L < it, small 
fluctuations will tend to drive the value of 9 towards ir/2. 

5.3.2 The sphere 

For the sphere distribution, we recall that the two allowed orientations of the string lead 
to quite different behaviors and that a confining potential appears in the 9o = n/2 case. 
The results of the stability analysis for the two orientations are as follows. 

• #o = 0. In this case, we have the Schrodinger potentials 



V x (u;u ) 



V e (u;u ) 



8u 8 - 18m 6 + 11m 4 + [2u 2 (u 2 - 1) - l]u 2 + u 2 (u 2 - 1) 
4u 2 {u 2 - I) 2 

2u 6 -u 4 + [Qu 2 (u 2 - 1) - l]u 2 - 3u 2 (u 2 - 1) 
4u 2 {u 2 - l) 2 ' 



(5.21) 



V y (u;u ) 



8u 8 - 18m 6 - [8u 2 (u 2 - 1) - llK + [Qu 2 (u 2 - 1) - l}u 2 - 3u 2 (u 2 - 1) 

4u 2 {u 2 - l) 2 



30 



and the value of the endpoint xq reads 

x = Vl - 2k' 2 K(k') , (5.22) 

where k' is the complementary modulus defined in (4.29). We have the general behavior 
(3.12) and 

7T 37T 

x (u ) ~ - - —(uo - 1) + O(u - I f . (5.23) 

Since V x and Vg are positive for all values of the parameter u$, the solution is stable under 
transverse and angular perturbations. On the other hand, V y starts from a negative value 
and turns positive. Repeating the analysis of section 5.1, we find that a longitudinal 
instability occurs for uq below the the critical value uq c given in the first of (4.31). 

• 9 = n/2. Now, the Schrodinger potentials are given by 

2u e -u 4 -uj 
V x (u;u ) = s , 

XL 

V e (u;u ) = -1, (5.24) 

2m 6 - m 4 - 2uiu 2 + 3ui 
V y {u;u ) = — . 

Hi 

and the value of the endpoint xq reads 



'9P — 1 

*o=^^— K(Jfc), (5.25) 

where k is the modulus defined in (4.34). We have the general behavior (3.12) and 

x (uo) ~ In + o ((uq - 1) ln(«o - 1)) • (5.26) 

Since 14 and are positive for all values of u , the solution is stable against transverse 
and longitudinal perturbations, again in accordance with the general results of section 
3.3. On the other hand, Vg has a constant negative value which in particular implies that 
the infinite- well approximation is exact. Examining the behavior of the lowest eigenvalue 
u)q, we find that it is an increasing function starting at negative values for u = and 
changing sign at a critical value uo c with 

u 0c ~ 1.14 , L c ~ 1.7 , E c ~ 0.22 . (5.27) 

That is, when the quark- ant iquark separation becomes larger than L c , small fluctuations 
in 6 destabilize the classical solutions and the resulting potential should not be trusted. 
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Since the configurations giving rise to a linear potential correspond to separations larger 
than L c , the confining behavior is spurious and instead we have a screened Coulomb 
potential with screening length given by the second of (5.27). Moreover, we may show 
that for u < u 0c more states become negative. Indeed, setting u>n( u o) = an d using 
(5.26) we find that the value uo c>n in which the n-th energy eigenvalue becomes zero is 

u 0c>n -1 + 16 e ^ 2 " +1 ) , (5.28) 

valid, practically, for all n ^ 1. 

The upshot of this analysis is that there is no stable confining branch and that both 
potentials are of the screened Coulomb type, with comparable screening lengths. 

5.4 Special points 

We have seen that, once we cross the critical value for u from above, the lowest eigenvalue 
of the fluctuations in the y-direction becomes negative and the only way for it to turn back 
to positive values is the appearance of another extremum of the length at a different value 
of Uq. However, as we have already mentioned the singularity structure of the fluctuation 
equations (3.6) changes when Uq = u m [ n . This isolated point corresponds to zero length 
and energy and to two straight strings stuck together. It is easily seen that we can have a 
positive-definite spectrum of fluctuations and therefore perturbative stability. However, 
for fluctuations with a parameter u infinitesimally larger than u min , the spectrum has a 
single negative eigenvalue. As we shall see, the apparent paradox is resolved by the fact 
that perturbation theory breaks down when applied to points in the vicinity of u = u m i n . 
This shows that the special points with u = u min are really of measure zero in all physical 
processes, in the sense that no conclusion reached at these points remains approximately 
correct when we move, even infinitesimally, away from them. 

This can be easily argued for the case of the longitudinal fluctuations for the non-extremal 
D3-brane. When u = 1 we see from (5.4) that the potential becomes 

u 8 - 1 

V y (u; 1) = 2^— , u > 1 , (5.29) 

Hi 

whereas from (5.6) we see that the Schrodinger equation is defined in the entire positive 
half-line. Hence the spectrum of fluctuations is positive and any perturbative analysis 
around uq — 1 necessarily breaks down. 
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A similar argument holds for the case of angular fluctuations for the disc and the tra- 
jectory corresponding to 9q = 0. In that case the Schrodinger problem can be expressed 
explicitly in terms of the new variable x in (3.8) as 

1 

u = — — , < x < oo . 5.30 
smha; 

The Schrodinger potential in terms of the variable x is given by 

V «« 0) = k{ l -Jk) ■ (531) 

It can be shown that the solution is given in terms of hypergeometric functions and that 
the spectrum is continuous with a gap, i.e. oj 2 > 1/4. Since we know from the analysis 
above that the spectrum is actually negative below the value given approximately in 
(5.16), we conclude, as before, that the perturbative analysis around u = necessarily 
breaks down. 

Next, we demonstrate explicitly all details of this phenomenon in the particular case 
of the Coulomb branch for the sphere and for the trajectory with 9q = 0, in which 
case longitudinal fluctuations are unstable. If u = 1, the Schrodinger potential of the 
longitudinal fluctuations can be expressed explicitly as a function of the variable x of 
(3.8) which, for our case, reads 

u = — *— , ^ x ^ - . (5.32) 

SIM 2 



it is given by 



V y (x; 1) = -4- - — ^- - 7 , (5.33) 
y sm 2 x 4 cos 2 a; 4 



and falls into the class of Poschl-Teller potentials of type I. The corresponding differential 
equation has a complete set of orthogonal solutions given by 

tf n (ar) = V4n + 5 sin 2 x cos 1/2 x P^ /m (cos 2x) , n = 0, 1, . . . , (5.34) 

where are the Jacobi polynomials of n-th order. The respective eigenvalues are 

u 2 n = An 2 + lOn + 6 , n = 0,l,2,.... (5.35) 

This is a positive-definite spectrum, showing that small fluctuations do not destabilize 
this special point at which u — 1, L — and E = 0. Consider next a small deviation 
from the value u = 1. We will show that using (3.32) to compute the correction to 
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the potential leads to divergent integrals for the corrections to the energy eigenvalues in 
(5.35). For this, we need the expressions 



du 
dV„ 



= 3 



du 
dx 



du 



u =l 



UQ = 1 



u =l 



sin 2 x 



- sin 2 x I tan 4 x , 
2 ' 



1 



2 sinx 



;8 + tan 4 a;) , 



(5.36) 



3 

"47 ' 



where we have used (5.21) and (5.23) (equivalent to computing the integral in (3.35) with 
Uq — 1). From these we may explicitly compute the right hand side of (3.32), which is of 
the form 

1 



5V = (uq — l)v(x) , v(x) 



+ 0((x-n/2)-*) , 



(5.37) 



(x - tt/2) 6 

where v(x) is a complicated function, with the indicated singular behavior. From the 
expression (5.34) for the solutions to the unperturbed problem we find that, near x = n/2, 
^ ~ (tt/2 — x) 1 / 2 . Hence the integral f^ 2 dx ^ 2 5V diverges, thus proving the breakdown 
of perturbation theory near u — 1. 



6 Discussion 



In this paper we have examined the perturbative stability of string configurations dual 
to flux tubes between static quark- ant iquark pairs in J\f = 4 SYM at finite temperature 
and at the Coulomb branch. The motivation for our study was the fact that the quark- 
antiquark potentials computed via the AdS/CFT prescription for the above cases exhibit 
behaviors that are inconsistent with our field-theory expectations, namely (i) multiple 
branches of the potential, (ii) a heavily orient at ion- dependent screening length, and (iii) 
a linear confining behavior. Our stability analysis resolves the discrepancy by showing 
that the configurations corresponding to the upper branches of the potential are unstable 
against longitudinal perturbations while those giving rise to an orientation-dependent 
screening length and to a confining behavior are unstable under angular perturbations. 

The methods developed here can be extended to the more involved situation of string 
configurations in a boosted and/or rotating non-extremal D3-brane background, with the 
boost corresponding to a thermal medium moving with respect to the pair and the ro- 
tation corresponding to R-charge chemical potentials in the gauge theory. In such cases, 
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the metrics are no longer diagonal and the various fluctuations are not guaranteed to de- 
couple so that the general discussions of sections 2 and 3 must be modified. Nevertheless, 
we believe that an analytic treatment in these cases is also possible. 

Another potential application of these methods refers to Wilson-loop calculations in less 
super symmetric backgrounds. In fact, there are several examples [22, 23] where calcula- 
tions of the heavy quark-antiquark potential in backgrounds with Af = 1 supersymmetry 
yield, at large separations, a linear confining behavior which, in contrast to the J\f = 4 
case, is actually expected on physical grounds. Therefore, it would be particularly inter- 
esting to investigate whether the string configurations giving rise to such a behavior are 
stable. 
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A An analog from classical mechanics 

It is not often appreciated that the problem of calculating Wilson loops in the super- 
gravity approach, especially in cases where multiple branches of the solution appear, has 
striking similarities to a textbook problem in classical mechanics, namely that of deter- 
mining the shape of a thin soap film stretched between two rings (Plateau's problem). 3 
The main similarity of the two problems lies in the fact that, although the solution of 
the equations of motion is straightforward, the boundary conditions allow for multiple 
solutions and introduce a phase structure. Since a lot of insight for our problem can be 
gained by looking at this simpler situation, in what follows we give a modern pedagog- 
ical review of this mechanical analog. Details can be found in standard textbooks on 
variational methods (e.g. [25]), while a partial stability analysis has been done in [26]. 
3 See, however, [24] for a discussion of this analogy in the context of Wilson-loop correlators. 
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We consider a thin soap film stretched between two coaxial circular rings of unit radius, 
separated by a distance L. Neglecting gravity, we write the action as 



S = 



dt J do\doi^ 



-(xl + xl + xl)-l 



(A.l) 



where we have taken the mass density and the surface tension equal to 1 and \ respec- 
tively Here, (xi, X2, £3) are Cartesian coordinates on the embedding space, while (01, 02) 
and 7 Q/ 3 are the coordinates and the induced metric on the surface respectively. For a 
static, axially-symmetric configuration we introduce cylindrical coordinates (r, 0, z) in 
the embedding space, we use reparametrization invariance to set (a±, 02) = (z, <fi) and we 
choose the embedding r = r(z). Then the action reduces to 

L/2 



S = -2tt / dt / dzrVl + r' 2 . (A.2) 

J-L/2 



Independence of the Lagrangian from z leads to the first integral 



r 



= u , (A.3) 



where Uq is the value of r at the point where r'(z) = which by symmetry occurs at z — 0. 
Integrating (A.3) and imposing r'(0) = 0, we obtain the solution r(z) = u cosh(z/u ), 
first found by Euler and defining a surface of revolution known as the catenoid. The 
integration constant u is specified by the boundary condition r(±-|) = 1, which gives 

L = 2uq cosh -1 — . (A.4) 

The potential energy of the solution is 

E = txuqL (l + j- sinh = 2tt ^1 - u 2 + u 2 cosh -1 . (A.5) 

Note that these equations require ^ u ^ 1 . Besides the catenoid solution just 
discussed, there also exists the so-called Goldschmidt solution which describes two dis- 
connected circular films on the two rings, 4 with energy E = 2n. 

To examine the properties of these solutions, we note that since (A.4) gives L(0) = 
L{\) = and L"(uq) = ~ M0 ( 1 _^)3/2 < 0, L(uq) is a concave function of Uq with a single 
maximum for ^ uq ^ 1. Setting L'{uq) = 0, we obtain the transcendental equation 



1 - u 2 cosh" 1 — = 1 , (A.6) 



4 The existence of this solution is most clearly seen by using the parametrization (<ti, 172) = (r, 4>) and 
z = z(r), for which the first-order equation reads rz'/y/l + z' 2 = const, which is solved by z' = 0. 
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Figure 6: The three equilibrium shapes of the soap film, namely the shallow catenoid, 
the deep catenoid and the Goldschmidt solution, plotted at a natural scale for L = 1. 

whose solution is u 0c ~ 0.552 leading to the maximal length L c ~ 1.325. We must then 
distinguish between two cases. For L > L c , (A. 4) has no solution for u and the catenoid 
solution does not exist at all, leaving the Goldschmidt solution as the only available one. 5 
For L < L c , on the other hand, (A. 4) has two solutions for u , the largest (smallest) 
of which corresponds to a "shallow" ("deep") catenoid. 6 Meanwhile, the Goldschmidt 
solution exists as well, giving a total of three available solutions, shown in Fig. 6. To 
examine which one is energetically favored, we need to compare the energy of the shallow 
catenoid, E(uq) with u® > uq c , with that of the Goldschmidt solution, E = 2n. The 
former is a decreasing function of u , becoming equal to 2n at u = u 0c ~ 0.826 where 
the separation is L c ~ 1.055, and hence the lowest-energy solution is the shallow catenoid 
for uq > Uq c and the Goldschmidt solution for u < u 0c . To summarize, for L > L c the 
only possible solution is the Goldschmidt solution, while for L < L c all three solutions 
are available with the shallow catenoid being favored for L < L c and the Goldschmidt 
solution being favored for L > L c . This behavior is shown in Fig. 7. 

By now, the analogy with the Wilson-loop calculations in the main part of the paper 
should be obvious. Namely, the quantities uq, L and E correspond to the quantities 
denoted by the same symbols in the Wilson- loop context, the shallow catenoid, deep 
catenoid and the Goldschmidt solution correspond to the short string, the long string 
and the unbound configuration respectively and the critical values L c and L c of the 
separation correspond to the maximal and screening lengths respectively. Fig. 7(b) in 
particular is qualitatively similar to Figs. 1(b) and 4(b) in the main part of the paper. 

5 That is, when the two rings supporting the soap film are stretched a distance larger than L c apart, 
the soap film will break to form two flat circular films over each ring. 

6 For example, for L = -|, these solutions have u ~ 0.967 (almost cylindrical film) and u ~ 0.076 
(highly curved film) respectively. 
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Figure 7: Plots of L(u ) and E(L), with the solid dark, solid gray and dashed gray lines 
corresponding to stable, metastable and unstable catenoid solutions. In the plot of E(L), 
the dotted line corresponds to the Goldschmidt solution. 

To examine the stability of the catenoid solutions, we need to consider small fluc- 
tuations about these equilibrium surfaces. In [26], the stability analysis was carried 
out for perturbations normal to the surface. Here, we consider the most general type 
of perturbation, which we parametrize as follows. Writing the equilibrium surface as 
F(r,z) = r — r c \(z) = 0, where r c \(z) = u cosh(z / u ) is the classical solution, the unit 
normal to the surface is n = j^fp while the unit tangent vectors along and perpendicular 
to the azimuthal direction are and £ = n x 0, respectively. Explicitly, we have the 
following transformation between the two orthonormal frames with basis vectors (n, £, 4>) 
and (r, respectively 

~ z + r' cl f 



n 



r — r cl z 



Expressing the most general perturbation as 



(A.7) 



5f = 5nn + r5(p(p + <5££ = 5rf + r5(p(p + 5zz , 

we have, in terms of our original variables, 

6£(t, <7i, cr 2 ) - r' cl (ai)5n(t, a u a 2 ) 



z = a 1 J r 



^cl(o"i) + 



Sn(t, (Ji, cr 2 ) + r^(<7i)^(i, o x , a 2 ) 



>/!+*) 



= cr 2 + <7i, <7 2 ) , 



(A.8) 



(A.9) 



From this form we can see that this is an 5*0(2) transformation with rotation angle 
related to the variable 01 as cos# = 1/ cosh(^/« ). Since this is local in space and global 
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in time the kinetic energy remains diagonal in this new basis. Substituting in the original 
action (A.l) and expanding in powers of the perturbations, we find that the zeroth-order 
term gives the classical action, the first-order term vanishes by the equations of motion 
and periodicity of S(f> in a 2 , and the second-order term reads 

S 2 = -±-J dt J ' d^da^ulid^nf + (d 2 5n) 2 } 

_ < cos h 2 ^ (sh 2 + 5i 2 + ul cosh 2 

2 Uq \ Uq 



, Sn 5f + sinh — 5n5f 

cosh 2 ^ V 2 u 

UQ 



cosh £1 



dUSnSn + sinh — ( -d^i 2 + cosh —5£d 2 c 
u V 2 u 



+ u 2 Q cosh — (d^d^cf) - d^d^cj)) \ . (A. 10) 

u J 

Although the various perturbations appear to be coupled, the calculation of the equations 
of motion reveals that they actually decouple due to an extensive cancellation of terms. 
For the tangential perturbations, the equations of motion are just 

8'£ = 5'<j> = 0. (A.ll) 

For the normal perturbations, we rename (a±, a 2 ) — > (z, 0), we define u = ^, we separate 
variables according to 

Sn(t, u, 0) = <5>(u)e~ int e im * , (A. 12) 

and we end up with the Sturm-Liouville equation (3.7) with 

2 

p = 1 , r = g m2 ? Q — cosh 2 u , uo = u Q , (A. 13) 

cosh u 

subject to the following boundary conditions at the endpoints 

K^HH 11 " 1 ^ ' (A14) 

where we have used (A. 4). To investigate the stability of the catenoid solutions, we want 
to determine the sign of u 2 for the lowest-energy solution to the differential equation (3.7) 
with (A. 13) in terms of Uq. Although the study of the evolution of uo 2 as a function of Mo is 
in general a hard task, we can obtain useful information by considering the corresponding 
zero-mode problem i.e. solving (3.7) with uo 2 = 0. In this case, the transformation x = 
tanhw turns (3.7) into an associated Legendre equation with the general solution given 
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by a linear combination of P 1 m (tanhw) and Q™(tanhu), m = 0, ±1, ±2, .... However, 
the boundary conditions (A. 14) further restrict us to m = 0, for which only the solution 
proportional to Qi(tanhu) is acceptable. Therefore, the zero- mode solution reads 

$o(«) = iVQxOnhw) = jV(l - wtanhw) , (A.15) 

where N is a normalization constant. Imposing the boundary condition (A. 14) leads then 
to the transcendental equation (A. 6), that is, to the condition determining the point u 0c 
at which the length L{uq) attains its maximum! Therefore, the equation of motion of the 
normal perturbations has a zero mode if and only if Uq attains the value Woe that marks 
the boundary between the deep and shallow catenoid solutions. This result tremendously 
simplifies our problem, as it implies that if we determine the evolution of the lowest 
eigenvalue as we move infinitesimally from uq c , we will in fact have determined its sign 
in the whole regions to the left and to the right of uq c . Again, one may appreciate the 
similarity with the case of longitudinal perturbations in the main part of the paper. 

To carry out this investigation, it is convenient to transform the Sturm-Liouville problem 
to a Schrodinger one by applying the transformation (3.8) which, for our case, reads 

x = sinhu, V(x) = (1 + x 2 ) 1/A <$>{u) . (A.16) 

This way, we obtain the Schrodinger equation (3.9) where the potential, for general values 
of m, is given by 

•r t / \ m 2 X 2 +6 /A i rr\ 

= —i ~ K^Tf • < A - 17 > 



and the boundary conditions at the endpoints read 



2 



*(±x ) = , x = x {u ) = °- . (A.18) 

The potential is depicted in Fig. 8(a) and clearly does not support any negative energy 
states for \m\ ^ 2, which is consistent with the discussion above. For m = it is deep 
enough to support one negative energy state if its size is restricted to a finite interval 
x G [— xo,xq], with x < x 0c . For \m\ = 1 there is no negative or zero energy state 
respecting the boundary conditions as we have seen, even though the potential has a 
negative part. In the new variables, the normalized zero- mode eigenfunction occurring 
for m = is 

«-<'> - 1"s -L"^ 374 ' 1 + * 2)1/4 i 1 - vmt sinh " ') • (A - 19 » 
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(a) (b) (c) 

Figure 8: (a) Schrodinger potentials for the normal perturbations with angular quantum 
numbers m = (solid), m = ±1 (dashed) and m = ±2 (dotted), (b) Evolution of the 
lowest eigenvalue u 2 with uq, calculated using the shooting method and Mathematica's 
[27] NDSolve routine. The slope of the curve at u = u 0c is found to be 4.497, in perfect 
agreement with (A. 21). (c) Plots of the normalized ground-state wavefunction ^o(x) for 
uq = 0.8 > uo c (solid) and uq = 0.3 < woe (dashed). 

To determine the change in uj 2 as u deviates from uq c , we may use perturbation theory. 
Noting that the potential (A. 17) is wo^independent while the endpoints ±xq given in 
(A. 18) are ito-dependent, we find that the appropriate perturbation-theory formula is 

5uj 2 = 2x , (u 0c )* 2 (x 0c )V(x 0c ) 5u , (A.20) 

as is read off from (3.33), appropriately adapted to a Schrodinger problem for a real 
wavefunction in the range [— xo c ,^oc] and with the boundary condition in (A. 18). In our 
case xq c ~ 1.509. To calculate 5u 2 , we insert (A. 17) and (A. 19) into (A.20), making 
repeated use of (A. 6) to simplify the resulting expressions. When the smoke clears out, 
we find 

Suj 2 = l2u °\ 5u ~ 4.497 5u . (A.21) 
3-5< 

Hence, as we move to the right (left) of u 0c , the eigenvalue u 2 becomes positive (nega- 
tive) and therefore the shallow catenoid is stable while the deep catenoid is unstable, as 
expected. This behavior is confirmed by a numerical analysis, shown in Fig. 8. 
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